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Abstract 

One of the objectives of designing feature selection learning algorithms is to obtain classifiers that depend on a small number 
of attributes and have verifiable future performance guarantees. There are few, if any, approaches that successfully address the two 
goals simultaneously. Performance guarantees become crucial for tasks such as microarray data analysis due to very small sample 
sizes resulting in limited empirical evaluation. To the best of our knowledge, such algorithms that give theoretical bounds on the 
future performance have not been proposed so far in the context of the classification of gene expression data. In this work, we 
investigate the premise of learning a conjunction (or disjunction) of decision stumps in Occam's Razor, Sample Compression, and 
PAC-Bayes learning settings for identifying a small subset of attributes that can be used to perform reliable classification tasks. 
We apply the proposed approaches for gene identification from DNA microarray data and compare our results to those of well 
known successful approaches proposed for the task. We show that our algorithm not only finds hypotheses with much smaller 
number of genes while giving competitive classification accuracy but also have tight risk guarantees on future performance unlike 
other approaches. The proposed approaches are general and extensible in terms of both designing novel algorithms and application 
to other domains. 
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Index Terms 

Microarray data classification. Risk bounds. Feature selection. Gene identification. 

I. Introduction 

An important challenge in the problem of classification of high-dimensional data is to design a learning algorithm that can 
construct an accurate classifier that depends on the smallest possible number of attributes. Further, it is often desired that 
there be realizable guarantees associated with the future performance of such feature selection approaches. With the recent 
explosion in various technologies generating huge amounts of measurements, the problem of obtaining learning algorithms 
with performance guarantees has acquired a renewed interest. 

Consider the case of biological domain where the advent of microarray technologies JEisen and Brownlll999llLipshutz et al 



199911 have revolutionized the outlook on the investigation and analysis of genetic diseases. In parallel, on the classification 



front, many interesting results have appeared aiming to distinguish between two or more types of cells, (e.g. diseased vs. 
normal, o r cells with different types of cancers) based on gene expression data in the case of DNA microarrays (see, for 
instance, JAIon et aI.L 1199911 for results on Colon Cancer, JGoIub et all 1199911 for Leukaemia). Focusing on very few genes to 
give insight into the class association for a microarray sample is quite important owing to a variety of reasons. For instance, a 
small subset of genes is easier to analyze as opposed to the set of genes output by the DNA microarray chips. It also makes 
it relatively easier to deduce biological relationships among them as well as study their interactions. An approach able to 
identify a very few number of genes can facilitate customization of chips and validation experiments- making the utilization 
of microarray technology cheaper, affordable, and effective. 

In the view of a diseased versus a normal sample, these genes can be considered as indicators of the disease's cause. 
Subsequent validation study focused on these genes, their behavior, and their interactions, can lead to better understanding of 
th e disease. Some att empts in this direction have yielded interesting results. See, for instance, a recent algorithm proposed 
by IWang et al.l 1200711 involving the identification of a gene subs et based on impor t ance r anking and subsequently combinations 



of genes for classification. Another example is the approach of 'Tibshir ani et al.l 11200311 based on nearest s hrunken centroids 



Some kernel based approaches such as the BAHSIC algorithm [Song et al.L 1200711 and their extensions (e.g., JShah and Corbeil , 
2010(1 for short time-series domains) have also appeared. 



The tradi tional method s used for classifying high-dimensional da ta are often characterized as either "filters" (e.g. IFurev et al 



|2000, Wan g et al.L 1200711 or "wrappers" (e.g. JGuvon et aI.L 1200211 ) depending on whether the attribute selection is performed 



independent of, or in conjunction with, the base learning algorithm. 
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Despite the acceptable empirical results achieved by such approaches, there is no theoretical justification of their performance 
nor do they come with a guarantee on how well will they perform in the future. What is really needed is a learning algorithm 
that has pwvably good performance guarantees in the presence of many irrelevant attributes. This is the focus of the work 
presented here. 

A. Contributions 

The main contributions of this work come in the form of formulation of feature selection strategies within well established 
learning settings resulting in learning algorithms that combine the tasks of feature selection and discriminative learning. 
Consequently, we obtain feature selection algorithms for classification with tight realizable guarantees on their generalization 
error The proposed approaches are a step towards more general learning strategies that combine feature selection with the 
classification algorithm and have tight realizable guarantees. We apply the approaches to the task of classifying microarray 
data where the attributes of the data sample correspond to the expression level measurements of various genes. In fact the 
choice of decision stumps as learning bias has in part motivated by this application. The framework is general and extensible 
in a variety of ways. For instance, the learning strategies proposed in this work can readily be extended to other similar tasks 
that can benefit from this learning bias. An immediate example would be classifying data from other microarray technologies 
such as in the case of Chromatin Immunoprecipitation experiments. Similarly, learning biases other than the conjunctions of 
decision stumps, can also be explored in the same frameworks leading to novel learning algorithms. 

B. Motivation 

For learning the class of conjunctions of features, we draw motivation from the guarantee that exists for this class in the 
following form: if there exists a conjunction, that depends on r o ut of the n inpu t attributes and that correctly classifies a 



training set of m examples, then the greedy covering algorithm o nHausslen kl988il will find a conjunction of at most r In 



iWL 



attributes that makes no training errors. Note the absence of dependence on the number n of input attributes. The method is 
guaranteed to find at most r In m attributes and, hence, depends on the number of available samples m but not on the number 
of attributes n to be analyzed. 

We propose learning algorithms for building small conjunctions of decision stumps. We examine three approaches to obtain 
an optimal classifier based on this premise that mainly vary in the coding strategies for the threshold of each decision stump. 
The first two approaches attempt to do this by encoding the threshold either with message strings (Occam's Razor) or by 
using training examples (Sample Compression). The third strategy (PAC-Bayes) attempts to examine if an optimal classifier 
can be obtained by trading off the sparsitjl^ of the classifier with the magnitude of the separating margin of each decision 
stump. In each case, we derive an upper bound on the generalization error of the classifier and subsequently use it to guide the 
respective algorithm. Finally, we present empirical results on the microarray data classification tasks and compare our results 
to the state-of-the-art approaches proposed for the task including the Sup port Vector Machine (SVM ) coupled with feature 



selectors, and Adaboost. The preliminary results of this work appeared in IMarchand and Shahl 1200511 . 



C. Organization 

Section |ll] gives the basic definitions and notions of the learning setting that we utilize and also characterizes the hypothesis 
class of conjunctions of decision stumps. All subsequent learning algorithms are proposed to learn this hypothesis class. 
Section Hill proposes an Occam's Razor approach to learn conjunctions of decision stumps leading to an upper bound on the 
generalization error in this framework. Section |IV] then proposes an alternate encoding strategy for the message strings using 
the Sample Compression framework and gives a corresponding risk bound. In Section |V] we propose a PAC-Bayes approach 
to learn conjunction of decision stumps that enables the learning algorithm to perform an explicit non-trivial margin-sparsity 
trade-off to obtain more general classifiers. Section |VI] then proposes algorithms to learn in the three learning settings proposed 
in Sections [Till |IV] and [V] along with a time complexity analysis. Note that the learning (optimization) strategies proposed in 
Section |Vl] do not affect the respective theoretical guarantees of the learning settings. The algorithms are evaluated empirically 
on real world microarray datasets in Section Ivnl Section IVIIII presents a discussion on the results and also provides an analysis 
of the biological relevance of the selected genes in the case of each dataset, and their agreement with published findings. Finally, 
we conclude in Section HXl 

II. Definitions 

The input space X consists of all ?i-dimensional vectors x = [xi, . . . , x„) where each real-valued component Xi e [Ai, Bi] 
for i = 1, . . .n. Each attribute Xi for instance can refer to the expression level of gene i. Hence, Ai and Bi are, respectively, 
the a priori lower and upper bounds on values for Xi. The output space y is the set of classification labels that can be assigned 
to any input vector x E X. We focus here on binary classification problems. Thus 3^ = {0, 1}. Each example z ~ (x, y) is an 

' This refers to the number of decision stumps used. 



input vector x € X with its classification label y E y chosen i.i.d. from an unknown distribution D on X x y. The true risk 
R{f) of any classifier / is defined as the probability that it misclassifies an example drawn according to D: 

R{f) 1^' P^^,y)^D (/(x) ^y)^ E(,,,)^,,/(/(x) ^ y) 

where I {a) = 1 if predicate a is true and otherwise. Given a training set S = {zi, . . . ,Zm} of m examples, the empirical 
risk Rs (/) on S, of any classifier /, is defined according to: 

_j rn 
' i=l 

The goal of any learning algorithm is to find the classifier with minimal true risk based on measuring empirical risk (and other 
properties) on the training sample S. 

We focus on learning algorithms that construct a conjunction of decision stumps from a training set. Each decision stump 
is just a threshold classifier defined on a single attribute (component) x^. More formally, a decision stump is identified by an 
attribute index k G {1, . . . , n}, a threshold value t £ M, and a direction d G { — 1, +1} (that specifies whether class 1 is on 
the largest or smallest values of Xk)- Given any input example x, the output r'^j^i^) of a decision stump is defined as: 

fc / N dcf J 1 if {xk — t)d > 
rtdW - I if ixk-t)d<0 

We use a vector k ^ (fci, . . . , fc|k|) of attribute indices kj 6 {1, . . . ,n} such that fci < fc2 < ■ • • < fc|k| where |k| is the 
number of indices present in k (and thus the number of decision stumps in the conjunction) B Furthermore, We use a vector 
t = (ifcj , tfe2 , . . . , ifc|j^| ) of threshold values and a vector d = {dk^ ,dk2, ■ ■ ■ , dk^^^ ) of directions where kj G {1, . . . ,n} for 
j G {1, . . . , |k|}. On any input example x, the output C^j(x) of a conjunction of decision stumps is given by: 

^w .dcf f 1 if r^,,,(x) = l Vj-Gk 
^*'*^^^-\ if 3jGk:r^-^,^.(x)=0 

Finally, any algorithm that builds a conjunction can be used to build a disjunction just by exchanging the role of the positive 
and negative labeled examples. In order to keep our description simple, we describe here only the case of a conjunction. 
However, the case of disjunction follows symmetrically. 

III. An Occam's Razor Approach 

Our first approach towards learning the conjunction (or disjunction) of decision stumps is the Occam's Razor approach. 
Basically, we wish to obtain a hypothesis that can be coded using the least number of bits. We first propose an Occam's Razor 
risk bound which will ultimately guide the learning algorithm. 

In the case of zero-one loss, we can model the risk of the classifier as a binomial. Let Bin(K, m, r) be the the binomial tail 
associated with a classifier of (true) risk r. Then Bin(K, to, r) is the probability that this classifier makes at most n errors on 
a set of m examples: 

Bin {n, m, r) 1^' ^ h\ r\l - r)™"* 

The binomial tail inversion Bin (k, to,, 5) then gives the largest ri sk value that a classifier can have while still h aving a probability 
of at least 5 of observing at most n errors out of m examples IILangfordi l2005[ Blum and Langfordi 1200311 : 



Bin (k, 771, 5) = sup {r : Bin (k, m, r) > S} 

From this definition, it follows that Bin (77T,i?5(/), m, (5) is the smallest upper bound, which holds with probability at least 
1 — (5, on the true risk of any classifier / with an observed empirical risk Rs{f) on a test set of m examples: 

V/: Prs^D" (i?(/) <BS(TOi?s(/),™,<5)) >l-6 

Our starting point i s the Occam's razor bound of iLangfordl 1200511 which is a tighter version of the bound proposed 
bv lBlumer et alj il987ll . It is also more general in the sense that it applies to any prior distribution P over any countable class 



of classifiers. 

-Although it is possible to use up to two decision stumps on any attribute, we limit ourselves here to the case where each attiibute can be used for only 
one decision stump. 



Theorem 1 dLangfordl 1200511 ). For any prior distribution P over any countable class J- of classifiers, for any data-generating 



distribution D, and for any 5 S (0, 1], we have: 

Prs-D- |v/ e F: R{f) < BS(mi?5(/), m, P{f)5) \>l-5 



The proof (available in ILangford, 1200511 ) directly follows from a straightforward union bound argument and from the fact 
that '^f^jrP{f) ~ 1- To apply this bound for conjunctions of decision stumps we thus need to choose a suitable prior P for 
this class. Moreover, Theorem[T]is valid when J^fej^^if) — 1- Consequently, we will use a subprior P whose sum is < 1. 

In our case, decision-stumps' conjunctions are specified in terms of the discrete-valued vectors k and d and the continuous- 
valued vector t. We will see below that we will use a finite-precision bit string a to specify the set of threshold values t. Let 
us denote by P(k, d, a) the prior probability assigned to the conjunction C^^ described by (k, d, a). We choose a prior of 
the following form: 

P(k,d,cr) = -77rTP(|k|)-^5k,d(c^) 
l|k|j ^ 

where gk.dlc) is the prior probability assigned to string a given that we have chosen k and d. Let A1(k, d) be the 
set of all message strings that we can use given that we have chosen k and d. If I denotes the set of all 2" possi- 
ble attribute index vectors and X'k denotes the set of all 2^^^ binary direction vectors d of dimension |k|, we have that 
EkeiEdePkE<xe.M(k,d)^(k,d,a) < 1 whenever ELoM'^) < 1 ^nd E^eA4(k.d) 9k,d(^) < 1 Vk, d. 

The reasons motivating this choice for the prior are the following. The first two factors come from the belief that the final 
classifier, constructed from the group of attributes specified by k, should depend only on the number |k| of attributes in this 
group. If we have complete ignorance about the number of decision stumps the final classifier is likely to have, we should 
choose p{d) = l/{n + 1) for d £ {0, 1, . . . , n}. However, we should choose a p that decreases as we increase d if we have 
reasons to believe that the number of decision stumps of the final classifier will be much smaller than n. Since this is usually 
our case, we propose to use: 

p(|k|) = A(|k| + l)-2 
TT- 

The third factor of P(k, d, a) gives equal prior probabilities for each of the two possible values of direction dj. 

To specify the distribution of strings .gk.dCo"), consider the problem of coding a threshold value t € [a,b] C [^,5] where 
[A, B] is some predefined interval in which we are permitted to choose t and where [a, b] is an interval of "equally good" 
threshold values^ We propose the following diadic coding scheme for the identification of a threshold value that belongs to 
that interval. Let / be the number of bits that we use for the code. Then, a code of I bits specifies one value among the set A; 
of threshold values: 



A,'^=^^ 



27-1 
1- ■' 



2'^ 






We denote by Ai and Bi, the respective a priori minimum and maximum values that the attribute i can take. These values are 
obtained from the definition of data. Hence, for an attribute i S k, given an interval [a,, hi] C [Ai, Bi] of threshold values, we 
take the smallest number Z, of bits such that there exists a threshold value in A;, that falls in the interval [a^, hi]. In that way, 
we will need at most [log2((i?i — Ai)/{bi — ai))\ bits to obtain a threshold value that falls in [a^, hi]. 

Hence, to specify the threshold for each decision stump i e k, we need to specify the number li of bits and a /j-bit string 
Si that identifies one of the threshold values in Aj^. The risk bound does not depend on how we actually code a (for some 
receiver). It only depends on the a priori probabilities we assign to each possible realization of a. We choose the following 
distribution: 

.gk,d(CT) =^ ffk,d('i,si,...,/|k|,s|k|) (1) 

= l[C{h)-2-'^ (2) 

iek 

where: 

C(a) =' A(a + i)-2 VaeN (3) 

The sum over all the possible realizations of a gives 1 since Ei^i ^ ^ = tt^/G. Note that by giving equal a priori probability 
to each of the 2'' strings Si of length U, we give no preference to any threshold value in A/.. 

The distribution ( that we have chosen for each string length li has the advantage of decreasing slowly so that the risk 
bound does not deteriorate too rapidly as k increases. Other choices are clearly possible. However, note that the dominant 
contribution comes from the 2^'' term yielding a risk bound that depends linearly in li. 

By a "good" threshold value, we mean a threshold value for a decision stump that would cover many negative examples and very few positive examples 
(see the learning algorithm). 



With this choice of prior, we have the following theorem: 
Theorem 2. Given all our previous definitions and for any 5 G (0, 1], we have: 

V.^^^ (Vk,d,a: R[C^^) <WR[mRs{C^^lm^^-^^^^^0^y^ >l~6 

Finally, we emphasize that the risk bound of Theorem [2] used in conjunction with the distribution of messages given by 
gk,d{<^), provides a guide for choosing the optimal classifier. Note that the above risk bound suggests a non-trivial trade-off 
between the number of attributes and the length of the message string used to encode the classifier Indeed the risk bound 
may be smaller for a conjunction having a large number of attributes with small message strings (i.e., small /^s) than for a 
conjunction having a small number of attributes but with large message strings. 

IV. A Sample Compression Approach 



The basic idea of the Sample compression framework IKuzmin and Warmuthl 1200711 is to obtain learning algorithms with 



the property that the generated classifier (with respect to some training data) can often be reconstructed with a very small 
subset of training examples. More formally, a learning algorithm A is said to be a sample-compression algorithm iff there 
exists a compression function C and a reconstruction function TZ such that for any training sample S = {zi, . . . , z™} (where 
Zi = {xi,yi)), the classifier A{S) returned by A is given by: 

A{S) = TZ{C{S)) \/S e{X X y)"" 

For a training set S, the compression function C of learning algorithm A outputs a subset Zj of S, called the compression set, 
and an information message a, i.e., (zi, a) = C(zi, . . . , z™). The information message <t contains the additional information 
needed to reconstruct the classifier from the compression set Zi. Given a training sample S, we define the compression set Zi 
by a vector of indices i such that i = (ii, 12, ... , J|i|), with ij G {1, . . . , m}\/j and ii < 12 < . . . < i|i| and where |i| denotes 
the number of indices present in i. 

When given an arbitrary compression set Zi and an arbitrary information message a, the reconstruction function 7?. of a 
learning algorithm A must output a classifier The information message a is chosen from a set 7W(zi) that consists of all the 
distinct messages that can be attached to the compression set Zi. The existence of this reconstruction function TZ assures that 
the classifier returned by A{S) is always identified by a compression set Zj and an information message tr. 

In sample compression settings for learning decision stumps' conjunctions, the message string consists of the attributes and 
directions defined above. However, the thresholds are now specified by training examples. Hence, if we have |k| attributes 
where k is the set of thresholds, the compression set consists of |k| trainin g examples (one per thre s hold). 



Our starting point is the following generic Sample Compression bound llMarchand and Sokolovai 1200511 : 



Theorem 3. For any sample compression learning algorithm with a reconstruction function TZ that maps arbitrary subsets of 
a training set and information messages to classifiers: 

Ps~D"> {Vi e X, a e M{Zi) -. R{TZ{a, Z;)) < e(a, Z;, \i\)} > 1 - S 
e(cr,zi,|j|) = 1-cxp 



where 

-1 



"^-|i|-|Jl 



In ( , 1 + In 



ln( ^^-l+ln' ^ 



PMizoi<y)J VC(|i|)C(lJl)^ 



(4) 



and C is defined by Equation \3\ 



Now, we need to specify the distribution of messages (^x(Zi)(f)) for the conjunction of decision stumps. Note that in 
order to specify a conjunction of decision stumps, the compression set consists of one example per decision stump. For each 
decision stump we have one attribute and a corresponding threshold value determined by the numerical value that this attribute 
takes on the training example. 

The learner chooses an attribute whose threshold is identified by the associated training example. The set of these training 
examples form the compression set. Finally, the learner chooses a direction for each attribute. 

The subset of attributes that specifies the decision stumps in our compression set Zi is given by the vector k defined in the 

previous section. Moreover, since there is one decision stump corresponding to each example in the compression set, we have 

i| = |k|. Now, we assign equal probability to each possible set |k| of attributes (and hence thresholds) that can be selected 

from n attributes. Moreover, we assign equal probability over the direction that each decision stump can have (+1, —1). Hence, 

we get the following distribution of messages: 



Pa.(zo('^)=(|J^|) •2^'^fVa (5) 

Equation |5] along with the Sample Compression Theorem completes the bound for the conjunction of decision stumps. 

V. A PAC-Bayes Approach 

The Occam's Razor and Sample Compression, in a sense, aim at obtaining sparse classifiers with minimum number of stumps. 
This sparsity is enforced by selecting the classifiers with minimal encoding of the message strings and the compression set in 
respective cases. 

We now examine if by sacrificing this sparsity in terms of a larger separating margin around the decision boundary (yielding 
more con fidence) can l ead u s to classifiers with smaller generalization error. The learning algorithm is based on the PAC-Bayes 
approach JMcAllesten, 1199911 that aims at providing Probably Approximately Correct (PAC) guarantees to "Bayesian" learning 
algorithms specified in terms of a prior distribution P (before the observation of the data) and a data-dependent, posterior 
distribution Q over a space of classifiers. 

We formulate a learning algorithm that outputs a stochastic classifier, called the Gibbs Classifier Gq defined by a data- 
dependent posterior Q. Our classifier will be partly stochastic in the sense that we will formulate a posterior over the threshold 
values utihzed by the decision stumps while still retaining the deterministic nature for the selected attributes and directions for 
the decision stumps. 

Given an input example x, the Gibbs classifier first selects a classifier h according to the posterior distribution Q and then 
use h to assign the label h{x.) to x. The risk of Gq is defined as the expected risk of classifiers drawn according to Q: 

RiGg) "^^ I]h^QR{h) = Eh^QEi^^y-^^oHH^) / y) 



Our starting point is the PAC-Bayes theorem JMcAllesten, l2003i iLangfordl l2005l ISeegen, 1200211 that provides a bound on the 
risk of the Gibbs classifier: 

Theorem 4. Given any space Ti of classifiers. For any data-independent prior distribution P over %, we have: 

V.^^ (VQ : kl(i?,(Gg)||i?(Gg)) < ^l^mEll}^^ > 1 _ , 
where KL(Q||P) is the Kullback-Leibler divergence between distribution^] Q and P: 

and where ]d{q\\p) is the Kullback-Leibler divergence between the Bernoulli distributions with probabilities of success q and 
p: _ 

kl(g|b)1^^1n^ + (l-g)ln-^ 
P 1 -P 

This bound for the risk of Gibbs classifiers can easily be turned into a bound for the risk of Bayes classifiers Bq over the 
posterior Q. Bq basically performs a majority vote (under measure Q) of binary classifiers in Ti. When Bq misclassifies an 
example x, at least half of the binary classifiers (under measure Q) misclassifies x. It follows that the error rate of Gq is at 
least half of the eiTor rate of Bq. Hence R{Bq) < 2R{Gq). 

In our case, we have seen that decision stump conjunctions are specified in terms of a mixture of discrete parameters k and 
d and continuous parameters t. If we denote by Pk,d(t) the probability density function associated with a prior P over the 
class of decision stump conjunctions, we consider here priors of the form: 



As before, we have that: 



whenever X]e=oP(^) = ^■ 

The factors relating to the discrete components k and d have the same rationale as in the case of the Occam's Razor 
approach. However, in the case of the threshold for each decision stumps, we now consider an explicitly continuous uniform 

■^Here Q{h) denotes the probability density function associated with Q, evaluated at h. 
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prior. As in the Occam's Razor case, we assume each attribute value Xk to be constrained, a priori, in [Ak , Bk] such that Ak 
and Bk are obtained from the definition of the data. Hence, we have chosen a uniform prior probability density on [Ak, Bk] 
for each tk such that fc e k. This explains the last factors of Pk,d(t). 

Given a training set S, the learner will choose an attribute group k and a direction vector d deterministically. We pose the 
problem of choosing the threshold in a similar manner as in the case of Occam's Razor approach of Section |lll] with the only 
difference that the learner identifies the interval and selects a threshold stochastically. For each attribute Xk G [Ak , Bk] : fc e k, 
a margin interval [ak,bk] ^ [Ak,Bk] is chosen by the learner A deterministic decision stump conjunction classifier is then 
specified by choosing the thresholds values tk G [ak, bk] uniformly. It is tempting at this point to choose tk = {ak+bk)/2 Vfc £ k 
(i.e., in the middle of each interval). However, the PAC-Bayes theorem offers a better guarantee for another type of deterministic 
classifier as we see below. 

Hence, the Gibbs classifier is defined with a posterior distribution Q having all its weight on the same k and d as chosen 
by the learner but where each tk is uniformly chosen in [ak,bk]- The KL divergence between this posterior Q and the prior 
P is then given by: 



k|/ P(|k|)y ^ \bk-ak 

In this limit when [ak, bk] = [Ak, Bk] V/c G k, it can be seen that the KL divergence between the "continuous components" 
of Q and P vanishes. Furthermore, the KL divergence between the "discrete components" of Q and P is small for small values 
of |k| (whenever p(|k|) is not too small). Hence, this KL divergence between our choices for Q and P exhibits a tradeoff 
between margins (bk — ak) and sparsity (small value of |k|j /or Gibbs classifiers. Theorem |4] suggests that the Gq with the 
smallest guarantee of risk R{Gq) should minimize a non trivial combination of KL{Q[[P) and Rs{Gq). 

The posterior Q is identified by an attribute group vector k, a direction vector d, and intervals [ak, bk] Vfc G k. We refine the 
notation for our Gibbs classifier Gq to reflect this. Hence, we use G^jJ where a and b are the vectors formed by the unions 
of flfeS and 6feS respectively. We can obtain a closed-form expression for ^sCG^b) ^y ^^^^ considering the risk P(x,y)(G'^b) 
on a single example (x, y) since Rs{G^^) = E(x,i;)~sP(x,y)(G^fJ). From our definition for Q, we find that: 



^(x,,)(G^C) = (l-2y) 



IKnM-y 



.fcGk 



(6) 



where: 



<^Ux) 



dcf 



if (.T < a and d = +1) or (b < x and d = —1) 
\f a < X <b and d = +1 
if a < X < 6 and d = — 1 

1 if (6 < a; and d = +1) or [x < a and d ~ —1) 



x — a 
b—a 
b—x 
b—a 



Note that the expression for P(x,y)(C'td) i^ identical to the expression for R{x.y}{G^^) except that the piece-wise linear 
functions a J' j^ (xk) are replaced by the indicator functions I{{xk — tk)dk > 0). 

The PAC-Bayes theorem provides a risk bound for the Gibbs classifier G^jJ. Since the Bayes classifier B^^ just performs 
a majority vote under the same posterior distribution as the one used by G^jJ, it follows that: 

Note that B^^ has an hyperbolic decision surface. Consequently, B^^ is not representable as a conjunction of decision 
stumps. There is, however, no computational difficulty at obtaining the output of i?^{J(x) for any x G A". We now state our 
main theorem: 

Theorem 5. Given all our previous definitions, for any S G (0, 1], and for any p satisfying X]e=oP(^) ~ ^' ^^ have, with 
probability atleast 1 — 6 over random draws of S ^ D"^: 

(Vk,d,a,b: R{G^t)<s^^p{^- kl(i?s(G^^)lk) < ^} ) 
where 

Furthermore: i?(B^{J) < 2R{G^^) Vk, d, a, b. 
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VI. The Learning Algorithms 



Having proposed the theoretical frameworks attempting to obtain the optimal classifiers based on various optimization criteria, 
we now detail the learning algorithms for these approaches. Ideally, we would like to find a conjunction of decision stumps 
that minimizes the respective risk bounds for each approach. Unfortunately, this cannot be don e efficiently in all cases since 
this pr oblem is at least as hard as the (NP-hard) minimum set cover problem as mentioned by iMarchand and Shawe-Tavlor 
1 200211 . Hence, we use a set covering greedy heuristic. It consists of choosing the decision stump i with the largest utility Uf'-' 
where: 

uf'' = m~pm (8) 

where Qi is the set of negative examples covered (classified as 0) by feature i, Ri is the set of positive examples misclassified 
by this feature, and p is a learning parameter that gives a penalty p for each misclassified positive example. Once the feature 
with the largest Ui is found, we remove Qi and Ri from the training set S and then repeat (on the remaining examples) until 
either no mo re negative examples are present or t hat a maximum number s of features has been reached. This heuristic was 
also used by IMarchand and Shawe-Tavloij l2002ll in the context of a sample compression classifier called the set covering 
machine. For our sample compression approach (SC), we use the above utility function Ui'^ . 

However, for the Occam's Razor and the PAC-Bayes approaches, we need utility functions that can incorporate the opti- 
mization aspects suggested by these approaches. 



A. The Occam's Razor learning algorithm 

We propose the following learning strategy for Occam's Razor learning of conjunctions of decision stumps. For a fixed U 
and ?7, let N be the set of negative examples and P be the set of positive examples. We start with N' — N and P' = P. Let 
Qi be the subset of N' covered by decision stump i, let Ri be the subset of P' covered by decision stump i, and let k be 
the number of bits used to code the threshold of decision stump i. We choose the decision stump i that maximizes the utility 
jjOccam defined as: 



u: 



v -k 



Occam dgf \Qt\ _ |-^»l 

N' ^ P 

where p is the penalty suffered by covering (and hence, misclassifying) a positive example and ?] is the cost of using li bits 
for decision stump i. Once we have found a decision stump maximizing Ui, we update N' — N' ~ Qi and P' = P' — Ri and 
repeat to find the next decision stump until either A^' = or the maximum number v of decision stumps has been reached 
(early stopping the greedy). The best values for the learning parameters p, 77, and v are determined by cross-validation. 



B. The PAC-Bayes Learning Algorithm 

Theorem |5] suggests that the learner should try to find the Bayes classifier B^^ that uses a small number of attributes {i.e., 
a small |k|), each with a large separating margin {hk — ak), while keeping the empirical Gibbs risk Rs{G^^) at a low value. 
As discussed earlier, we utilize the greedy set covering heuristic for learning. 

In our case, however, we need to keep the Gibbs risk on S low instead of the risk of a deterministic classifier Since the 
Gibbs risk is a "soft measure" that uses the piece-wise linear functions erf ^ instead of the "hard" indicator functions, we 
cannot make use of the hard utility function of Equation [8] Instead, we need a "softer" version of this utility function to take 
into account covering (and erring on) an example partly. That is, a negative example that falls in the linear region of a erf ^ is 
in fact partly covered and vice versa for the positive example. 

Following this observation, let k' be the vector of indices of the attributes that we have used so far for the construction of 
the classifier. Let us first define the covering value C{G^^) of G^^^ by the "amount" of negative examples assigned to class 
by G^^^: 



C{G^^^) ''^^ E (1-y) 






We also define the positive-side error £{G^y^) of G^^ as the "amount" of positive examples assigned to class : 



(x,y)es 



jek' 



We now want to add another decision stump on another attribute, call it i, to obtain a new vector k" containing this new 
attribute in addition to those present in k'. Hence, we now introduce the covering contribution of decision stump i as: 
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and the positive-side error contribution of decision stump i as: 

Typically, the covering contribution of decision stump i should increase its "utility" and its positive-side error should 
decrease it. Moreover, we want to decrease the "utility" of decision stump i by an amount which would become large 
whenever it has a small separating margin. Our expression for KL{Q\\P) suggests that this amount should be proportional 
to \n{{Bi — Ai)/{bi — ai)). Furthermore we should compare this margin term with the fraction of the remaining negative 
examples that decision stump i has covered (instead of the absolute amount of negative examples covered). Hence the covering 
contribution C^-^{i) of decision stump i should be divided by the amount M^■^^ of negative examples that remains to be 
covered before considering decision stump i: 

(x,y)eS jek' 

which is simply the amount of negative examples that have been assigned to class 1 by G\^. If P denotes the set of positive 
examples, we define the utility U^^{i) of adding decision stump i to G^J* as: 



Ae^'i ^ \p\ 



rk'd/-\ dcf l-ab \^) „^ab 



^ab 



where parameter p represents the penalty of misclassifying a positive example and ?/ is another parameter that controls the 
importance of having a large margin. These learning parameters can be chosen by cross-validation. For fixed values of these 
parameters, the "soft greedy" algorithm simply consists of adding, to the current Gibbs classifier, a decision stump with 
maximum added utility until either the maximum number v of decision stumps has been reached or all the negative examples 
have been (totally) covered. It is understood that, during this soft greedy algorithm, we can remove an example (x, y) from S 
whenever it is totally covered. This occurs whenever Oigk' '^a bi^i) ~ ^■ 

Hence, we use the above utility function for the PAC-Bayes learning strategy. Note that, in the case of U^^ and jjOccam^ 
we normalize the number of covered and erred examples so as to increase their sensitivity to the respective rj terms. 

1) Time Complexity Analysis: Let us analyze the time complexity of this algorithm for fixed p and rj. For each attribute, 
we first sort the m examples with respect to their values for the attribute under consideration. This takes 0{m\ogm) time. 
Then, we examine each potential ai value (defined by the values of that attribute on the examples). Corresponding to each ai, 
we examine all the potential bi values (all the values greater than a,)- This gives us a time complexity of 0{m?). Now if k is 
the largest number of examples falling into [a^, hi], calculating the covering and error contributions and then finding the best 
interval [ai,fei] takes 0{km^) time. Moreover, we allow k G 0{m) giving us a time complexity of 0{'m?) for each attribute. 
Finally, we do this over all the attributes. Hence, the overall time complexity of the algorithm is 0{nm^). Note, however, 
that for microarray data, we have n » m (hence, we can consider m?' to be a constant). Moreover once the best stump is 
found, we remove the examples covered by this stump from the training set and repeat the algorithm. Now, we know that 
greedy algorithms of this kind have the following guarantee: if there exist r decision stumps that covers all the m examples, 
the greedy algorithm will find at most r\x\{m) decision stumps. Since we almost always have r G 0(1), the running time of 
the whole algorithm will almost always be G 0{nriv' log(?Ti)). The good news is, since n » m, the time complexity of our 
algorithm is roughly linear in n. 

2) Fixed-Margin Heuristic: In order to show why we prefer a uniformly distributed threshold as opposed to the one fixed 
at the middle of the interval [ai,6i] for each stump i, we use an alternate algorithm that we call the fixed margin heuristic. 
The algorithm is similar to the one described above but with an additional parameter 7. This parameter decides a fixed 
margin boundary around the threshold, i.e. 7 decides the length of the interval [ai, hi]. The algorithm still chooses the attribute 
vector k, the direction vector d and the vectors a and b. However, the a^'s and 6,'s for each stump i are chosen such that, 
\bi — flil = 27. The threshold ti is then fixed in the middle of this interval, that is U = ^"'"^ '' . Hence, for each stump i, the 
interval [ai, hi] = [U — 7, ti + 7]. For fixed p and 7, a similar analysis as in the previous subsection yields a time complexity 
of 0{nm? log(TO)) for this algorithm. 
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Data Set 


SVM 


SVM+ 


^s 


SVM+rfe 


Adaboost | 


Name 


ex 


Genes 


Krrs 


Errs 


S 


Errs 


S 


Itrs 


Errs 


Colon 


62 


2000 


12.8±1.4 


14.4±3.5 


256 


15.4±48 


128 


20 


15.2±2.1 


B MD 


34 


7129 


I3.2±l 


7.2±2.6 


32 


10.4±2.4 


64 


20 


9.8±1.1 


C MD 


60 


7129 


28.2±2.2 


23.1±2.8 


1024 


28.2±2.2 


7129 


50 


21.2±2.4 


Leuk 


72 


7129 


21.3±1.4 


14±2.8 


64 


21±3.2 


256 


20 


17.8±1.8 


Lung 


52 


918 


8.8±1.3 


6.8±1.9 


64 


7.2±1.8 


32 


1 


2.4±1.4 


BreastER 


49 


7129 


15.3±2.4 


10.3±2.7 


256 


11.2±2.8 


256 


50 


9.8±1.7 



TABLE I 
Results of SVM, SVM coupled with Golub's feature selection algorithm (filter), SVM with Recursive Feature Elimination 

(WRAPPER) AND ADABOOST ALGORITHMS ON GENE EXPRESSION DATASETS. 



Data Set 


Occam 


SC 1 


Name 


ex 


Genes 


Errs 


S 


bits 


Errs 


S 


Colon 
B MD 
C_MD 
Leuk 
Lung 
BreastER 


62 
34 
60 

72 
52 
49 


2000 
7129 
7129 
7129 
918 
7129 


23.6±1.2 
17.2±1.8 
28.6±1.8 
27.8±1.7 
21.7±1.1 
25.4±1.2 


1.8±.6 
1.2±.8 
2.6±1.1 
2.2±.8 
1.8±1.2 
3.2±.6 


6 

3 
4 
6 

5 
2 


18.2±1.8 
17.2±1.3 
29.2±1.1 
27.3±1.7 
18±1.3 
21.2±1.5 


1.2±.6 
1.4±.8 
1.2±.6 
1.4±.7 
1.2±.5 
1.4±.5 



TABLE II 

Results of the proposed Occam's Razor and Sample Compression learning algorithms on Gene Expression datasets. 



VII. Empirical Results 

The proposed approaches f or learning conjun ctions of decisio n stumps were teste d on the six real-world binary microarray 
datas ets viz. the colon tumo r JAlon et aL ^ 199911. the Leuka emia IIGolub et al.ll 199911 . the B_MD and C MP Medulloblastomas 
data JPomeroy et al.L l2002ll the Lun2 llOarber et al.L I2OOII1 . and the BreastER data llWest et al.[ I2OOII1 . 

The colon tumor data set llAlon et al.L 1 199911 provides the expression levels of 40 tumor and 22 normal colon tissues measured 
for 6500 hum an genes. We use th e set of 2000 genes identified to have the highest minimal intensity across the 62 tissues. The 
Leuk data set JGolub et al.l 1199911 provides the expression levels of 7129 human genes for 47 samples of patients with Acute 
Lymphob lastic Leukemia (ALL ) and 25 samples of patients with Acute Myeloid Leukemia (AML). The B_MD and C_MD 
data sets OPomeroy et all 1200211 are microarray samples containing the expression levels of 7129 human genes. Data set B_MD 
contains 25 classic and 9 desmoplastic medulloblastomas whereas data set C_MD contains 39 medulloblastomas survivors and 
21 treatment failures (non-survivors). The Lung dataset consis t s of g ene expression levels of 918 genes of 52 patients with 39 
Adenocarcinoma and 13 Squamous Cell Cancer iGarber et all 1200 ill . This data has some missing values which were replaced 
by zeros. Finally, the BreastER dataset is the Breast Tumor data of IWest et al.l 1200111 used with Estrogen Receptor status 
to label the various samples. The data consists of expression levels of 7129 genes of 49 patients with 25 positive Estrogen 
Receptor samples and 24 negative Estrogen Receptor samples. 

The number of examples and the number of genes in each data are given in the "ex" and "Genes" columns respectively under 
the "Data Set" tab in each table. The algorithms are referred to as "Occam" (Occam's Razor), "SC" (Sample Compression) and 
"R\C-Bayes" (PAC-Bayes) in Tables HIl to Iv] They utilize the respective theoretical frameworks discussed in Sections HiH IIVI 
and [V] along with the respective learning strategies of Section |VT] 

We have compared our learning algorithm with a linear-kernel soft-ma rgin SVM tr a ined b oth on all the attributes (gene 
expressions) and on a subset of attributes chosen by the filter method of iGolub et al.l 11199911 . The filter method consists of 
ranking the attributes as function of the difference between the positive-example mean and the negat ive-example mean a nd then 
use only t he first I attributes. T he resulting learning algorithm, named SVM+gs is the one used by iFurev et al.l 12000(1 for the 
same t ask. Guvon et al. 120021 claimed o btaining better results with the recursive feature elimination method but, as pointed 
out by lAmbroise and McLachlanI ll2002ll . their work contained a methodological flaw. We use the SVM recursive feature 
elimination algorithm with this bias removed and present these results as well for comparison (referred to as "SVM-i-rfe" in 
Table |I]i. Finally, we also compare our results with the state-o f-the-art Adaboost algorithm. For this, we use the implementation 
in the Weka data mining software IWitten and Franki 1200511 . 

Each algorithm was tested over 20 random permutations of the datasets, with the 5-fold cross validation (CV) method. Each 
of the five training sets and testing sets was the same for all algorithms. The learning parameters of all algorithms and the 
gene subsets (for "SVM-i-gs" and "SVM-i-rfe") were chosen from the training sets only. This was done by performing a second 
(nested) 5-fold CV on each training set. 

For the gene subset selecti on procedure of SV M-i-gs, we have considered the first I = 2' genes (for i = 0,1,...,12) ranked 
according to the criterion of iGolub et al.l 11 199911 and have chosen the i value that gave the smallest 5-fold CV error on the 
training set. The "Errs" column under each algorithm in Tables H] to |lll] refer to the average (nested) 5-fold cross-validation 
error of the respective algorithm with one standard deviation two-sided confidence interval. The "bits" column in Table HI] 
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Data Set 


PAC-Bayes | 


Name 


ex 


Genes 


S 


G-errs 


B-errs 


Colon 


62 


2000 


1.53±.28 


14.68±1.8 


14.65±1.8 


B MD 


34 


7129 


1.2±.25 


8.89±1.65 


8.6±1.4 


C MD 


60 


7129 


3.4±1.8 


23.8±1.7 


22.9±1.65 


Leuk 


72 


7129 


3.2±1.4 


24.4±1.5 


23.6±1.6 


Lung 


52 


918 


1.2±.3 


4.4±.6 


4.2±.8 


BreastER 


49 


7129 


2.6±1.1 


12.8±.8 


12.4±.78 



TABLE III 

Results of the PAC-Bayes learning algorithm on Gene Expression datasets. 



refer to the number of bits used for the Occam' Razor approach. The "G-errs" and the "B-errs" columns in Table Hill refer to 
the average nested 5-fold CV error of the optimal Gibbs classifier and the corresponding Bayes classifier with one standard 
deviation two-sided interval respectively. 

For Adaboost, 10, 20, 50, 100, 200, 500, 1000 and 2000 iterations for each datasets were tried and the reported results 
correspond to the best obtained 5-fold CV error. The size values reported here (the "S" columns for "SVMH-gs"and "SVM-nrfe", 
and "Itr" column for "AdaBoost" in Table |I]i correspond to the number of attributes (genes) selected most frequently by the 
respective algorithms over all the permutation runs|j Choosing, by cross-validation, the number of boosting iteration is somewhat 
inconsistent with Adaboost's goal of minimizing the empirical exponential risk. Indeed, to comply with Adaboost's goal, we 
should choose a large-enough number of boosting rounds that assures the convergence of the empirical exponential risk to its 
minimum value. However, as shown by IZhang and Yul ll2005ll . Boosting is known to overfit when the number of attributes 
exceeds the number of examples. This happens in the case of microarray experiments frequently where the number of genes 
far exceeds the number of samples, and is also the case in the datasets mentioned above. Early stopping is the recommended 
approach in such cases and hence we have followed the method described above to obtained the best number of boosting 
iterations. 

Further, Table |IV] gives the result for a single run of the deterministic algorithm using the fixed-margin heuristic described 
above. Table |V] gives the results for the PAC-Bayes bound values for the results obtained for a single run of the PAC-Bayes 
algorithm on the respective microarray data sets. Recall that the PAC-Bayes bound provides a uniform upper bound on the 
risk of the Gibbs classifier. The column labels refer to the same quantities as above although the errors reported are over a 
single nested 5-fold CV run. The "Ratio" column of Table IV] refers to the average value of {bk — ak)/{Bk — Ak) obtained 
over the decision stumps used by the classifiers over 5 testing folds and the "Bound" columns of Tables ITV] and [V] refer to the 
average risk bound of Theorem |5] multiplied by the total number of examples in respective data sets. Note, again, that these 
results are on a single permutation of the datasets and are presented just to illustrate the practicality of the risk bound and the 
rationale of not choosing the fixed-margin heuristic over the current learning strategy. 



A. A Note on the Risk Bound 

Note that the risk bounds are quite effective and their relevance should not be misconstrued by observing the results in just 
the current scenario. One of the most limiting factor in the current analysis is the unavailability of microarray data with larger 
number of examples. As the number of examples increase, the risk bound of Theorem 5 gives tighter guarantees. Consider, 
for instance, if the datasets for the Lung and Colon Cancer had 500 examples. A classifier with the same performance over 
500 examples (i.e. with the same classification accuracy and number of features as currently) would have a bound of about 12 
and 30 percent error instead of current 34.6 and 54.6 percent respectively. This only illustrates how the bound can be more 
effective as a guarantee when used on datasets with more examples. Similarly, a dataset of 1000 examples for Breast Cancer 

''There were no close ties with classifiers with fewer genes. 



Data Set 


Stumps:PAC-K 


ayes(fixed margin) 


Name 


ex 


Genes 


Size 


Errors 


Bound 


Colon 


62 


2000 


1 


14 


34 


B MD 


34 


7129 


1 


7 


20 


C MD 


60 


7129 


3 


28 


48 


Leuk 


72 


7129 


2 


21 


46 


Lung 


52 


918 


2 


9 


29 


BreastER 


49 


7129 


3 


11 


31 



TABLE IV 
Results of the PAC-Bayes Approach with Fixed-Margin Heuristic on Gene Expression Datasets. 
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Data Set 


StumpsiPAC-Bayes | 


Name 


ex 


Genes 


Ratio 


Size 


G-errs 


B-errs 


Bound 


Colon 


62 


2000 


0.42 


1 


12 


11 


33 


B MD 


34 


7129 


0.10 


1 


7 


7 


20 


C MD 


60 


7129 


0.08 


5 


21 


20 


45 


Leuk 


72 


7129 


0.002 


3 


22 


21 


48 


Lung 


52 


918 


0.12 


1 


3 


3 


18 


BreastER 


49 


7129 


0.09 


2 


11 


11 


29 



TABLE V 
An illustration of the PAC-Bayes risk bound on a sample run of the PAC-Bayes algorithm. 



with a similar performance can have a bound of about 30 percent instead of current 63 percent. Hence, the current limitation 
in the practical application of the bound comes from limited data availability. As the number of examples increase, the bounds 
provides tighter guarantees and become more significant. 

VIII. Analysis 

The results clearly show that even though "Occam" and "SC" are able to find sparse classifiers (with very few genes), they 
are not able to obtain acceptable classification accuracies. One possible explanation is that these two approaches focus on 
the most succinct classifier with their respective criterion. The Sample compression approach tries to minimize the number of 
genes used but does not take into account the magnitude of the separating margin and hence compromises accuracy. On the 
other hand, the Occam's Razor approach tries to find a classifier that depends on margin only indirectly. Approaches based on 
sample compression as well as minimum description length have shown encouraging results in various domains. An alternate 
explanation for their suboptimal performance here can be seen in terms of extremely limited sample sizes. As a result, the 
gain in accuracy does not offset the cost of adding additional features in the conjunction. The PAC-Bayes approach seems to 
alleviate these problems by performing a significant margin-sparsity tradeoff. That is, the advantage of adding a new feature is 
seen in terms of a combination of the gain in both margin and the empirical risk. This can be compared to the strategy used 
by the regularization approaches. The classification accuracy of PAC-Bayes algorithm is competitive with the best performing 
classifier but has an added advantage, quite importantly, of using very few genes. 

For the PAC-Bayes approach, we expect the Bayes classifier to generally perform better than the Gibbs classifier. This is 
reflected to some extent in the empirical results for Colon, C_MD and Leukaemia datasets. However, there is no means to 
prove that this will always be the case. It should be noted that there exist several different utility functions that we can use 
for each of the proposed learning approaches. We have tried some of these and reported results only for the ones that were 
found to be the best (and discussed in the description of the corresponding learning algorithms). 

A noteworthy observation with regard to Adaboost is that the gene subset identified by this algorithm almost always include 
the ones found by the proposed PAC-Bayes approach for decision stumps. Most notably, the only gene Cyclin Dl, a well 
known marker for Cancer, found for the lung cancer dataset is the most discriminating factor and is commonly found by both 
approaches. In both cases, the size of the classifier is almost always restricted to 1. These observations not only give insights 
into the absolute peaks worth investigating but also experimentally validates the proposed approaches. 

Finally, many of the genes identified by the^«fl|j PAC-Bayes classifier include some prominent markers for the corresponding 
diseases as detailed below. 



A. Biological Relevance of the Selected Features 

Table [Vll details the genes identified by the final PAC-Bayes classifier learned over each dataset after the parameter selection 
phase. There are some prominent markers identified by the classifier. Some of the main genes identified by the PAC-Bayes 
approach are the ones identified by previous studies for each disease — giving confidence in the proposed approach. Some 
of the discovered genes in this case include Human monocyte-de rived neutrophil-activating protein (MONAP) mRNA in the 
case of Colon Cancer dataset and oestrogen receptor in the case of Breast Cancer data, D79205 _at-Ribosomal protein L39, 
D83542_at-Cadherin-1 5 and U29195_at-NPTX2 Neuronal pentraxin II in the case of Medulloblastomas datasets B_MD and 
C_MD. Other genes identified have biological relevance, for instance, the identific ation of Adip sin, LAF -4 and HOXIC with 



1 2003 1 



reg ard to ALL/AML by our algori thm is in agreement with that of the findings of iChow et al.l [2001], Hiwatari et al 
and [Lawrence and LargmanI 11 199211 respectively and the studies that followed. 

Further, in the case of breast cancer. Estrogen receptors (ER) have shown to interact with BRCAl to regul ate VEGF 
transcription and secretion in breast cancer cells [Kawai et al., 2 002|l. These interacti ons are further investigated by iMa et al. 



II2OO5I1 . Further studies for ER have also been done. For instance. iMoggs et al.l 1200511 discovered 3 putative estrogen-response 
elements in Keratin6 (the second gene identified by the PAC-Bayes classifier in the case of BreastER data) in the context of 



This is the classifier learned after choosing the best parameters using nested 5-fold CV and trained on the full dataset. 



Gene(s) identified by PAC-Bayes Classifier 



Colon 
B_MD 
C MD 



Lung 
BreastER 



1. Hsa627 M26383-Human monocyte-deiived neutrophil-activating protein (MONAP) mRNA 
1. D79205_at-Ribosomal protein L39 

1. S7l824_at-Neuial Cell Adhesion Molecule, Phosphatidylinositol-Linked Isoform Piecursor 

2. D83542_at-Cadherin-15 

3. U29195_at-NPTX2 Neuronal pentraxin II 

4. X73358_s_at-HAES-1 mRNA 

5. L36069_at-High conductance inward rectifier potassium channel alpha subunit mRNA 

1. M84526_at-DF D component of complement (adipsin) 

2. U34360_at-Lymphoid nuclear protein (LAF-4) mRNA 

3. M16937_at-Homeo box cl protein, mRNA 

1. GENE221X-IMAGE_841641-cyclin Dl (PRADl -parathyroid adenomatosis 1) Hs82932 AA487486 

1. X03635_at,X03635- class C, 20 probes, 20 in all_X03635 5885 - 6402 
Human mRNA for oestrogen receptor 

2. L42611_f_at, L42611- class A, 20 probes, 20 in L42611 1374-1954, 
Homo sapiens keratin 6 isofonn K6e KRTHE mRNA, complete cds 



TABLE VI 
Genes Identified by the Final PAC-Bayes Classifier 



E2-responsive genes identified by microarray analysis of MDA-MD-23 1 cells that re-express ERq.. A n important role played 
by cytokeratins in cancer development is also widely known (see for instance iGusterson et al.l 11200511 '). 

Furthermore, the importance of MONAP in the case of colon c ancer and Adipsi n in the case of leukaemia data has further 
been confirmed by various rank based algorithms as detailed by ISu et alj ll2003ll in the implementation of "RankGene", a 
program that analyzes and ranks genes for the gene expression data using eight ranking criteria including Information Gain 
(IG), Gini Index (GI), Max Minority (MM), Sum Minority (SM), Twoing Rule (TR), t-statistic (TT), Sum of variances (SV) 
and one-dimensional Support Vector Machine (IS). In the case of Colon Cancer data, MONAP is identified as the top ranked 
gene by four of the eight criteria (IG, SV, TR, GI), second by one (SM), eighth by one (MM) and in top 50 by IS. Similarly, 
in the case of Leukaemia data, Adipsin is top ranked by IS, fifth by SM, seventh by IG, SV, TR, GI and MM and is in top 
50 by TT. These observations provides a strong vaUdation for our approaches. 

Cyclin as identified in the case of Lung Cancer datase t is a well known marker for cell d ivision whose perturbations are 
considered to be one of the major factors causing cancer JDriscoU et al.L 1 19991 iMasaki et al.l l2003il . 

Finally, the discovered genes in the case of Medulloblastomas are important with regard to the neuronal functioning (esp. 
S71824, U29195 and L36039) and can have relevance for nervous system related tumors. 



IX. Conclusion 

Learning from high-dimensional data such as that from DNA microarrays can be quite challenging especially when the aim 
is to identify only a few attributes that characterizes the differences between two classes of data. We investigated the premise of 
learning conjunctions of decision stumps and proposed three formulations based on different learning principles. We observed 
that the approaches that aim solely to optimize sparsity or the message code with regard to the classifier's empirical risk limits 
the algorithm in terms of its generalization performance, at least in the present case of small dataset sizes. By trading-off the 
sparsity of the classifier with the separating margin in addition to the empirical risk, the PAC-Bayes approach seem to alleviate 
this problem to a significant extent. This allows the R\C-Bayes algorithm to yield competitive classification performance while 
at the same time utilizing significantly fewer attributes. 

As opposed to the traditional feature selection methods, the proposed approaches are accompanied by a theoretical justification 
of the performance. Moreover, the proposed algorithms embed the feature selection as a part of the learning process itselfO 
Furthermore, the generalization error bounds are practical and can potentially guide the model (parameter) selection. When 
applied to classify DNA microarray data, the genes identified by the proposed approaches are found to be biologically significant 
as experimentally validated by various studies, an empirical justification that the approaches can successfully perform meaningful 
feature selection. Consequently, this represents a significant improvement in the direction of successful integration of machine 
learning approaches for use in high-throughput data to provide meaningful, theoretically justifiable, and reliable results. Such 
approaches that yield a compressed view in terms of a small number of biological markers can lead to a targeted and well 
focussed study of the issue of interest. For instance, the approach can be utilized in identifying gene subsets from the microarray 
experiments that should be further validated using focused RT-PCR techniques which are otherwise both costly and impractical 
to perform on the full set of genes. 

Finally, as mentioned previously, the approaches presented in this wor have a wider relevance, and can have significant 
implications in the direction of designing theoretically justified feature selection algorithms. These are one of the few approaches 
that combines the feature selection with the learning process and provide generalization guarantees over the resulting classifiers 
simultaneously. This property assumes even more significance in the wake of limited size of microarray datasets since it limits 
the amount of empirical evaluation that can be reliably performed otherwise. Most natural extensions of the approaches and 
the learning bias proposed here would be in other similar domains including other forms of microarray experiments such as 



^Note that lHuang and Chand I2007II proposed one such approach. However, they need multiple SVM learning runs. Hence, their method basically works 
as a wrapper 
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Chromatin Immunoprecipitation promoter arrays (chlP-Chip) and from Protein arrays. Within the same learning settings, other 
learning biases can also be explored such as classifiers represented by features or sets of features built on subsets of attributes. 
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